1 



Baryons on the lattice 

in 

Q : G.S. BalP 

o 



in 
o 



^Dept. of Physics & Astronomy, The University of Glasgow, Glasgow G12 8QQ, UK 



I comment on progress of lattice QCD techniques and calculations. Recent results on 
^ ' pentaquark masses as well as of the spectrum of excited baryons are summarizned and 
interpreted. The present state of calculations of quantities related to the nucleon structure 
and of electromagnetic transition form factors is surveyed. 
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O . In the past two years several narrow hadronic resonances have been discovered: new 
bottomonium and charmonium states, the and at least two new D* mesons. It is more 
than twenty years ago that new states with widths < 10 MeV have been seen last, in the T 
system. This is an exciting situation as only a rather small number of discovered hadrons, 
Q . including most notably quarkonia, the hydrogen of QCD, are narrow. What makes this 
^ ! new era of hadron spectroscopy even more interesting is that many of these previously 
^ ' overlooked states appear to require more QCD than is allowed within simple quark-model 
qq mesons or qqq baryons: at least some of the new states constitute the anti-thesis to 
^ ■ the "hydrogen of QCD". If confirmed in a high statistics experiment, the G"*" pentaquark 
. baryon and possibly other exotic baryons will add even more to this excitement. 

Lattice QCD is ideally positioned to compute the spectrum of reasonably stable hadrons 
as well as of non-perturbative properties relating to their internal structure. The above 
described experimental discoveries were paralleled by significant advances in lattice meth- 
ods, enabling computer simulations of QCD to become a precision predictive tool. In many 
areas there are also lessons to be learned from combining lattice studies with effective field 
theory (EFT) methods and/or model assumptions. 

We have also witnessed experimental progress in the study of the spectrum of excited 
baryons and of electromagnetic transition form factors. Generalised parton distributions 
(GPDs) of the proton are now being studied intensively. As a stable particle the nucleon 
lends itself to lattice studies. In the case of spin-independent structure functions it will 
be hard for lattice simulations to compete with the experimental precision. However, 
in addition to the theoretical satisfaction of verifying experimental measurements, this 
provides an ideal test ground for the methods and approximations employed in lattice 
studies. Once spin and transversity are included into the description of the nucleon, the 
experimental situation is far less clean and here there is real potential for lattice prediction 
rather than postdiction. 

I will describe the present state of the field, interpret recent pentaquark and baryon 
mass calculations and briefly survey studies on the form and structure of baryons. 
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Figure 1. The fraction of the nucleon spin car- 
ried by its quarks |2|. 
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Figure 2. The lowest two energy levels 
in the QQ system with sea quarks jH]- 



2. THE LATTICE: WHERE ARE WE? WHERE DO WE GO? 

QCD can be regularized by introducing a space-time lattice cut-off a. The QCD cou- 
pling and Uf quark masses whose values are not predicted by QCD should then be matched 
to reproduce nj + l experimental measurements of hadronic properties, for instance hadron 
masses. Everything else is a prediction and in this sense lattice QCD is a first principles 
approach. The confinement of colour implies that finite size effects are usually tiny, as 
long as the spatial box extent La ^ fn'^. We are fortunate to find that lattice spacings 
= 1-4 GeV are sufficiently small to allow for controlled continuum limit extrapo- 
lations, a ^ 0. This means that L <^ 100 is sufficient, which makes QCD tractable on 
computers. 

On a lattice with V = L^T points, the lattice Dirac operator is a huge matrix of 
dimension 12V^. The inversion of this operator represents the major computational task 
of lattice QCD and this makes simulations incorporating sea quarks expensive. The 
algorithmic cost explodes with small vr masses, oc l/(m^a)~'^. A smaller also requires 
a larger spatial lattice volume and the scaling behaviour, keeping m^^La fixed, is even 
worse: oc 1 / {mT^a)-'^ . This is the main reason why many simulations are performed in the 
quenched approximation including only valence quarks and neglecting the polarization of 
the QCD vacuum due to sea quarks. While this approximation violates unitarity and does 
not even qualify as a quantum field theory, light hadron masses seem to agree within 10 % 
with experiment PP, indicating that the main effect of quark loops can be absorbed into 
redefinitions of the bare parameters of the theory. The quenched approximation however 
goes terribly wrong at least in the scalar and pseudoscalar flavour-singlet sectors. 

In Figure H we display recent n/ = 2 QCD results obtained by the LHP and SESAM 
Collaborations [2] on the quark contribution AS to the proton spin, in the MS scheme 
at a scale /i = 2 GeV. The normalization is such that | = ^AS + Lg + J^, where Lq 
is the contribution from the quark angular momentum and Jg from the gluons. Similar 
results have been obtained by the QCDSF Collaboration [l]. In these simulations, m^^ > 
550 MeV. Obviously, for infinite quark masses we expect AE = 1. It is therefore not 
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surprising that the experimental value is overestimated and it is clear that smaller quark 
masses are absolutely essential to allow for a meaningful chiral extrapolation. Fortunately, 
with the advent of new Fermion formulations [3J that respect an exact lattice chiral 
symmetry, a reduction of the quark mass towards values ^ 180 MeV has become 
possible jH], albeit so far only in the quenched approximation. 

In many lattice calculations it is sufficient to calculate quark propagators that originate 
from a fixed source point. In these cases only one column of the inverse Dirac matrix needs 
to be calculated, naively reducing the effort by a factor V. In some cases diagrams with 
disconnected quark lines are needed. Examples are the physics of flavour singlet mesons, 
strong decays as well as parton distributions. In the latter case the complication can be 
avoided by assuming SU{2) isospin symmetry and only calculating differences between 
u and d quark distributions. Disconnected quark lines require all-to-all propagators and 
hence a complete inversion of the Dirac matrix appears necessary. This turns out to be 
prohibitively expensive in terms of memory and computer time. Fortunately, sophisticated 
noise reduced stochastic estimator techniques have been developed over the past few years 
and as a result tremendous progress was achieved. One such benchmark is the QCD string 
breaking problem, Q{r)Q(0) ^ B(r)B{0), where B = Qq and Q is a static quark. This 
represents one of the cleanest examples of a strong decay. Within both, the transition 
matrix element as well as the BB state all-to-all propagators are required. In Figure |21 we 
display the result of a recent SESAM Collaboration study |S| with Uf = 2, mq ^ mg and 
a ~ 0.085 fm. An extrapolation to physical light quark masses yields a string breaking 
distance Tc ~ 1.16 fm. The gap between the two states in the string breaking region is 
AE ^ 50 MeV and we are able to resolve this with a resolution of 10 standard deviations! 

In conclusion, all the long standing "killers" <^ rris, nj > and all-to-all propagators 
have been successfully tackled. However, we are still a few years away from overcoming 
combinations of two of these simultaneously and possibly up to a decade separates us from 
precision simulations of flavour singlet diagrams with realistically light sea quarks. Not 
all these ingredients are always required at the same time. While ever bigger computers 
are an absolute necessity, most of the recent progress would have been impossible without 
novel methods. The gain factor from faster computers was almost 5,000 over the past 15 
years. The factor from theoretical and algorithmic advances is harder to quantify. 

3. PENTAQUARKS 

QCD goes beyond the quark model and hence hadronic states that do not flt into a naive 
quark model of qq mesons and qqq baryons are of particular interest. Many of the observed 
hadrons will contain considerable higher Fock state components. Obviously, quantum 
numbers that are incomprehensible with a quark model meson or baryon interpretation 
provide us with the most clean-cut distinction. Such examples do exist in the Review of 
Particle Properties, namely the J^*" = 1 ^ mesons 7ri(1400) and 7ri(1600). The minimal 
conflguration required to obtain a vector state with positive charge either consists of two 
quarks and two antiquarks (tetraquark/molecule) or of quark, antiquark and a gluonic 
excitation (hybrid meson). These resonances are rather broad with widths P ~ 300 MeV. 
However, the ratio P/m is very much the same as for the established p(770) vector meson. 

Another clear indication of an rzquark > 3 nature would for instance be a baryonic state 
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with strangeness S* = +1. The minimal quark configuration in this case consists of five 
quarks (pentaquark): uudds. Over the past two years several experiments have presented 
evidence of a very narrow B"'"(1530) resonance [7,, with decay 6"*" — > pK^ and 0"*" — > nK^ . 
The parity has not yet been established. However the mass is about 100 MeV above the 
KN threshold and for = \ an S-wave decay is possible, which might be difficult to 
reconcile with a width F <^ 10 MeV. For |^ a P-wave is required, still a bit puzzling but 
less so. As the main decay channel does not require quark-antiquark pair creation one 
might hope to gain some insight from quenched lattice simulations and several attempts 
have been made ^8. 9 .10 11,12 13j . 

Two groups |9ll0j also investigated charmed pentaquarks and two studies |8lllj incor- 
porated the 1=1 sector, in addition to / = 0. Two groups |10|llj employed chiral 
overlap Fermions while the others used conventional Wilson-type lattice quarks. Only 
the Kentucky group lUJ went down to m^r ~ 180 MeV while all other vr masses were 
larger than 400 MeV. The Budapest- Wuppertal group [HI varied the lattice spacing and 
attempted a continuum-limit extrapolation. In all studies the negative parity mass came 
out to be lighter than the positive parity one, which is expected in the heavy quark limit. 

There are two crucial questions to be asked: what happens when realistically light quark 
masses are approached? Do we see resonant or scattering states? Resolving a resonance 
sitting on top of a tower of KN scattering states with different relative momenta appears 
rather hopeless at first. However, there are two discovery tools available: variation of the 
lattice volume and of the creation operator. By varying the volume (and the boundary 
conditions PU]) one will change the spectrum of KN scattering states as well as the 
coupling of a given operator to KN (the spectral weight). 

For the |^ state which can only decay into a P-wave, the mass of the scattering state 
will depend on the lattice size since the smallest possible non- vanishing lattice momentum 
is 7r/(aL). For ^ the volume dependence of the lowest scattering state mass will be weak, 
however, the scaling of the spectral weight with the volume provides us with additional 
information. 

It turns out that the situation on the lattice is at least as ambiguous as the one en- 
countered in experiment. To demonstrate this we display some Kentucky- Washington 
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results [TT] in Figure El It appears that the | state dominantly couples to an S'-wave 
KN . The I"*" displays the qualitative volume dependence of a P-wave, however, it does 
not share its non-interacting mass. It is most likely a P-wave scattering state. This 
interpretation is supported by the observed volume dependence of the spectral weight. 
At very small m^r the situation becomes further complicated by the fact that there is no 
axial anomaly in the quenched approximation. Hence the flavour singlet rj' is degenerate 
with the TT. As this contribution comes in with a negative spectral weight, it is sometimes 
labelled a "ghost". The I"*" state can contain such a KNr]' S'-wave (dashed lines). 

The ^ state becomes indistinguishable from a KN S'-wave as the quark mass is re- 
duced. This does not conclusively exclude the possible existence of a nearby resonant 
state which might only couple very weakly to the creation operator used in this particular 
study. In order to draw more definite conclusions a variation of the creation operator as 
well as of the volume appears necessary, which is a very ambitious project jT3]. If the 
pentaquark really was such a narrow resonance as some experiments suggest then maybe 
a lattice operator can be constructed that has a large overlap with this state but only a 
very small coupling to KN. 

Lattice studies of diquark interactions in a simplified, more controlled environment 
represent an alternative strategy to the brute force simulation of unstable states. A 
baryon with one static and two light quarks constitutes one such arena. One can of 
course also investigate multiquark interactions in the nonrelativistic limit of infinitely 
heavy quark masses. Such tetra- and penta-quark potentials have been studied recently 
by two groups |13|15|16j and the results should provide model builders with some insight. 
However, it is not clear how to relate these findings to the light quark limit in which chiral 
symmetry appears to play a bigger role than instantaneous confining forces. 

There exist quite a few narrow resonances very close to strong decay thresholds like 
the A(1405), the recently discovered X(3872) charmonium state and the ao//o(980) sys- 
tem. It is very conceivable that such states contain a sizable multiquark component. 
The question then arises if these are would-be quark model states or if these are true 
molecules/multiquark-states, that appear in addition to the quark model spectrum. A 
fantastic arena to address this was provided by the recently discovered (probably scalar) 
P'*(2317) and (probably axialvector) Z}*(2457) states. First lattice studies |17|18j have 
been performed, with somewhat contradictory interpretations of very compatible results. 
One might hope that a similar lattice effort will also be dedicated on the comparatively 
cleaner and easier question of tetraquarks as has been on pentaquarks. 

4. EXCITED BARYONS 

The spectrum of baryons has attracted renewed experimental and theoretical interest 
in recent years. There is the question if the states can be understood in terms of quark 
models and if so by what sort of interaction and assumptions. Do gluonic excitations or 
pentaquark components play a role for instance in the Roper resonance? What can we 
learn about quark-quark interactions within bound states? Quark model predictions are 
somewhat obscured as the corresponding decay widths set a limit on the precision that 
can be expected for the resulting masses. Are missing states really "missing" or are they 
just obliterated due to the presence of many very broad, overlapping resonances? Strongly 
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decaying hadrons also pose problems in lattice simulations. At present almost all calcu- 
lations of baryonic resonances have been performed within the quenched approximation 
in which these are stable and hence the problem is circumvented. 

This limitation can also be viewed as a virtue since most models suffer from an omission 
of quark pair creation effects too. Comparison with similar lattice results then allows to 
establish the validity range and applicability of a particular phenomenological model. 
One strength of lattice methods is that simulations are not limited to the quark mass 
parameters found in nature. Investigating the quark mass dependence of results is a 
powerful tool. In the limit of large quark masses one would expect to make contact with 
non-relativistic quark models while as m.^ overlap with chiral perturbation theory 
(xPT) predictions should be verifiable. 

Based on the assumptions that QCD bound states are mesons and baryons, that there 
is a mass gap and spontaneous chiral symmetry breaking at zero quark mass, an effective 
low energy chiral effective field theory (xEFT) can be derived in the spirit of the Born- 
Oppenheimer approximation. This will, to leading order, describe interactions between 
the (fast moving) massless Nambu-Goldstone pions and other hadrons. In nature quarks 
and thus pions are not massless and the leading mass corrections are formally of order 
m^/A^SB where A^^^ ^ 47rF^ > 1 GeV. The number of terms explodes at higher orders 
and predictive power is eventually lost, unless is sufficiently small to allow for an early 
truncation. 

Lattice simulations with sea quarks have so far been limited to unrealistically heavy 
pions, heavier than about 400 MeV. Only recently masses as low as 180 MeV have become 
possible jU] . To allow for a controlled extrapolation of lattice results to the physical region 
it is mandatory to establish an overlap between simulation data and xPT expectations. 
In general the size of this window will depend on the observable in question. Chiral 
lattice Fermion actions will make such a comparison with xPT cleaner. Only in this case 
an exact version of chiral symmetry can be formulated at finite lattice spacing. With 
other Fermion discretizations, strictly speaking, a comparison should only be attempted 
after extrapolating lattice results to the continuum limit. While even 400 MeV -C A^sb 
(modulo the ambiguity of vs. "<") such a pion is still doomed to "see" some of the 
internal structure of the proton, with an inverse charge radius of about 250 MeV. Hence 
it is doubtful that the quark and gluon nature of QCD can completely be ignored with 
such a "hard" pion probe. 

Naively, hadron masses are a polynomial in the quark mass, nig oc m^. However, pion 
loops give rise to a non-analytic functional dependence on the quark mass. For instance 
the nucleon mass is given by, 

mAr(m^) = m7v(0) + 02"^^ + a^'ml + 

with a renormalization scale A. The coefficients of the non-analytic terms can be related 
to phenomenological low energy constants. For instance 03 = —3g\/{32TTF^). Apart 
from such constants, 04 and a'^ contain terms oc m]^^ and a'^ a contribution oc 02. In the 
quenched approximation, the leading non-analytic term is not proportional to but to 
m^ln(m^/A), due to the rj' becoming an additional Goldstone pion. 

In Figure 0] we show a comparison [19l20j between Uf = 2 lattice data of the nucleon 
mass and the relativistic xPT expectation Eq. ^ [1^] . The lattice results were obtained 



e[{X) 04 + a^ln -7^ 
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Figure 4. Chiral extrapolation of the nucleon mass |19f2()j . 



by the CP-PACS and JLQCD Collaborations ^ as well as by the UKQCD and QCDSF 
Collaborations [22]; with (non-chiral) Wilson- type Fermions, on relatively fine lattices, 
a < 0.15 fm and large volumes, aLiriT^ > 5. The low energy parameters were fixed to 
phenomenological values and the fit comprises only of mi\f{0),C2 and e^(lGeV). The 
quantitative agreement between curve and data for > 600 MeV is accidental jlH] . A 
naive polynomial fit to the simulation data results in a nucleon mass much larger than the 
experimental value: lower order xPT only becomes applicable at smaller quark masses. 

EFTs are based on the separation of scales. If the xEFT however is regulated in 
dimensional regularization then vr loop-integrals can receive significant contributions from 
momenta q > A^sb- To enhance the convergence the Adelaide group |2S] suggested a 
"finite range regularization" approach which amounts to introducing a momentum cut-off 
which is then varied to achieve "model independence" . A hard cut-off can also be provided 
by lattice regularization of the chiral expansion |24|25j . Needless to say that all cut-off 
and scheme dependence will disappear at sufficiently high orders in the p-expansion. 

During the past four years we witnessed many lattice publications on the spectrum of 
excited baryons |2E]- While the extraction of a ground state mass is relatively straight 
forward, radial excitations either require high statistics and some confidence in the fitting 
procedures or the design of sophisticated, non-local creation operators ^7\. All but one 
study [2H] have been performed in the quenched approximation. Only the LHPC-UKQCD- 
QCDSF Collaborations [221, using the Wilson-clover action, attempted a continuum limit 
extrapolation. Other strategies were implementations of improved Wilson-type actions 
like the -D234 action by Lee et al. or the FLIC action by Melnitchouk et al. jSI] 
as well as a recent study with Wilson-clover Fermions jH21- The BGR Collaboration 
employed chirally improved Fermions [2Zj; the Riken-BNL group used chiral domain wall 
Fermions jSSj and the Taiwan jHlj as well as the Kentucky- Washington groups jO] made 
use of overlap Fermions. 

Early articles shared the observation of the positive parity state being much heavier than 
the Roper A^'(1440) resonance while the negative parity ground state was compatible with 
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the orbital excitation, A^*(1535). One explanation would be that the resonance observed 
in nature might have little overlap with the dominantly qqq state created on the lattice; 
pentaquark or gluon components might be necessary. Alternatively, maybe one should 
not take the exact position of a resonance with a width of 0(200 Me V) overly seriously. 
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Figure 5. The N , N' and A^* masses for Figure 6. Chirally improved Fermions |27j : 
overlap Fermions [HI at a ^ 1 GeV . 0.06 < m^/GeV^ < 0.65, a'^ ^ 1.33 GeV. 



New results obtained by the Kentucky-Washingon group [HI and the BGR Collabora- 
tion |2Z] at lighter pion masses, m.^- > 180 MeV and m.^. > 250 MeV, respectively, are 
however compatible with experiment. We display these in Figures and El Rather un- 
surprisingly the creation operator requires a node in its spatial wave function to produce 
a significant overlap with the radial excitation. This complication was circumvented in 
Ref. by a sophisticated fitting procedure while in Ref. |2Z1 such an adequate operator 
has been constructed, using a variational principle. Making contact with the light quark 
regime seems hopeless for > 500 MeV. However, the authors of two recent stud- 
ies |32|34j manage to extrapolate their results to the experimental values from such pion 
masses as well. Note that the BGR Collaboration only sees a clear signal of the radial 
excitation for m,r > 400 MeV while the Bayesian fitting procedure of the Kentucky group 
yields results at any quark mass. The Taiwan group ^34j in addition predicts the spectrum 
of doubly charmed baryons, with findings roughly compatible with earlier studies |35|36j 
as well as with the SELEX candidate(s). 

One might hope that in the near future the N' state can cleanly be disentangled from a 
P-wave rj'N scattering state or, at very small masses, from 5"- wave ttttN /rj'rj'N states. By 
studying the volume dependence and spectral weights the Kentucky group has taken steps 
in this direction. Most of the lattice studies include additional baryonic resonances where 
similar problems need to be addressed. In conclusion, we are close to an understanding of 
the transition between the heavy and light quark limits in the quenched approximation, 
an information invaluable to model builders. 
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5. FORM AND STRUCTURE 

Quite a few results on moments of GPDs, most notably from the QCDSF jHTj, SESAM 
and LHP Collaborations exist. A nice review of the state-of-the-art concerning spin- 
independent parton distributions and the axial charge can be found in Ref. [221 • The 
main problem here is an overestimation by about 60 % of {x)u-d if linearly extrapolated 
in m^, relative to experiment. It is not clear whether this difference will reduce as the 
quark mass is decreased or if this has to do with the non-chiral Fermion formulation used. 
This issue will be clarified in the near future [10] . 

There has also been progress in resolving the momentum dependence of electromagnetic 
7* A — >■ A transition form factors in a quenched study [H], in the region 0.1 GeV^ < < 
1.4 GeV^ The magnetic dipole form factor is significantly overestimated at large Q^, due 
to the unrealistically small charge radii of A and A. One might hope that such effects 
cancel in part from form factor ratios. Rem = Ge2/Gmi is fairly constant at -0.02(1), once 
extrapolated to the chiral limit. In contrast, Rcm = Gc2/Gmi decreases monotonously 
from -0.01(1) at 0.1 GeV^ down to -0.09(3) at > i QeV^. This behaviour is in good 
agreement with > 0.4 GeV^ CLAS data while it is hard to reconcile with the OOPS 
point RsM = [-6.1 ± 0.2 ± 0.5] % at ^ 0.13 GeVl 

6. SUMMARY 

Many lattice studies now include sea quarks. Within the quenched approximation, 
light quark masses close to the physical limit have been realised and the lattice provides 
a powerful tool for exploring the validity range of chiral expansions. Lattice pentaquark 
studies still yield ambiguous results. A systematic study using a large set of creation 
operators, of lattice volumes, spacings and quark masses is possible but ambitious. Latest 
lattice data suggest that the mass of the Roper resonance can be reproduced in the 
quenched approximation, thus indicating a non-exotic leading Fock component. A lot 
of progress has been made in understanding the structure of the nucleon and there is a 
strong push towards reducing the lattice quark masses, closer to the physical limit. 
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